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In this article we present all nonsingular upper triangular solutions to the constant quan- 
tum Yang-Baxter equation R^f 2 R l h lk ? R l ,? l l = R* 2 * 3 R.H 3 R{ ll l in three state case. The 
, upper triangular ansatz implies 729 equations for 45 variables. Fortunately many of the 

equations turned out to be simple allowing us to start breaking the problem into smaller 
ones. In the end we had a total of 552 solutions, but many of them were either inherited 
from two-state solutions or subcases of others. The final list contains 35 nontrivial solutions, 
most of them new. 
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1 Introduction 



In this paper we continue our work on systematically solving the constant quantum Yang-Baxter 
O '• equation 

Y\ ' jjk\k2 -nhkz -nhh 75^2^3 r>kil 3 jjhh f-i\ 

^ In our previous paper Q (for more details, see we solved completely the two-state problem, 

which involved 64 equations in 16 variables. 

In general (0) contains N G equations for iV 4 unknowns, so in the present case when N = 3 
we have 729 equations for 81 variables. At the moment this is too complicated for an exhaustive 
study and one must proceed step by step with suitable ansatze. The present ansatz is based on 
the two-state results: it turned out [||, ||, |j] that all two-state nonsingular solutions were either 

• Upper triangular: Rfj = 0, if j > I or j = I, i > k, or 

• Even weight: R$j = 0, if k + I ^ i + j (mod 2). 

(In thinking Rfj as an A^ 2 x A^ 2 matrix we use the convention that the right-hand indices 
define the block.) Thus upper triangularity turns out to be a good ansatz for finding two-state 
solutions, and we hope it is equally productive when the number of states is three. 
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2 Simple solutions 



Let us first of all note that the YBE does have some 'easy' classes of solutions. That is, the 
solutions are easy to construct from previous solutions or for some other reason. [This does not 
mean that the other related structures (e.g. algebras) are necessarily less interesting.] In this 
section we will describe two such classes and solutions belonging to them are not mentioned 
again when we discuss the results of our search. 

2.1 Solutions inherited from lower number of states 

A solution to (|l]) with a lower number of states can always be dressed to become a higher-state 
solution. There are at least two ways to do this. 



2.1.1 Diagonal dressing 

Let R be an M-state solution of ([l]) and N > M. Let \i be a selection of M numbers from 
{1, . . . , N}, and define the iV-state i?-matrix as follows: 



R 



kl 



f>kl 



when i, j, k,l G fi, 
otherwise. 



(2) 



When this is substituted into (|T|) we find further conditions on s as follows: 

R%j (•tmi'tmj SmkSml) — 0) 



(3) 



Vi, j, k, I G fi, m $l n, no sum. [It is interesting to note that the conditions are exactly like those 
obtained in a different context in M] A trivial solution to (||) is given by = 1, but there can 
also be others (depending on the form of R). 

In the present case we have M = 2 and N = 3 and there is only one index outside [i. As an 
example let us take \i = {1,3} and 



R 





■ ■ \ 


■ p 


1 — pq . 




Q 







(4) 



(for matrices we always write '.' in place of a '0', for better readability) and then the only 
remaining condition from (||) is S12S21 = S23S32. Thus we can extend (Q) to a three-sate solution 
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with the four additional parameters a, b, x, y: 

( 1 . • 



R 



P 



\ 



by 



1-pq 



ay 



(5) 
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2.1.2 Block dressing 

The starting point is as above but the new higher-sate solution is constructed as follows: 



kl 



R 



This ansatz leads to the conditions 



Rfj, when i,j,k,l £ fi, 
<5j fc Fj, when j,l G i,k £ pi 
G^5 l j, when i,k G fi, j,l £ \i 



(6) 



otherwise. 



{F®F)R 
(1 ® F)R(G <8> 1) 
R(G ® G) 
[F,G] 

The most important case of this is when G 
sufficient. An example of such an solution is 



R(F g) F), 

(G (g> 1)-R(1 <8> F), 

(G <g> G)R, 

0. 



(7) 



R := 



F 1 , in which case the first equation of (0) is 

\ 



1 . 1 




-1 . Cg 


. -1 . 




d e 


. 1 




Cg 










• h ■ 






96 








1 . -Cg 






. -1 






1 



(8) 
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which decomposes as 



R := 



V 



1 1 


-1 Cg 


. 1 


Cg 




1 -Cg 




1 



G :-- 



1 d 6 



G :-- 



96 dege 

96 



(9) 



2.2 Solutions generated by commuting matrices 

The following solutions can be said to be inherited from lower dimension, when dimension means 
the number of pairs of indices. (This hierarchial structure is more obvious when one considers 
extensions, see e.g. ((§.) 

The result is simply the following: Let {N(a),M(a)\a € 1} be a set of commuting N x N 
matrices, then it is easy to show that 



flg = £ N(a)fM(a) 



(10) 



is an iV-sate solution of (|T]). 

A commuting set of matrices can be simultaneously brought to the Jordan canonical form, 
but the Jordan form does not have to be diagonal. A two-state example is provided by 

/ 



R 



\ 



1 a 


b c \ 


. 1 


. b 




1 a 




■ l ) 



(11) 



it can be decomposed (this is not unique) as 

/ 1 • 



N(l) 



1 a 
. 1 



,M(1) 



1 



N(2) 



b c 
. b 



(12) 



When the number of states is three there is still more freedom and there are more varied solutions 
of this type, but we will not include such solutions in the list of results. An example is 

/ 1 . 1 . a 5 a 6 
1 . . . d 8 a 5 
1 



R := 



1 



Cg a s ag \ 
Cg dgdg 
Cg 



d 8 a 5 d$ciQ 
dja 5 

1 

1 

1 J 
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3 Symmetries 

It is well know that the set of equations (|l|) is invariant under the continuous transformation 

R -> k(Q ® Q)R(Q <g> Qy 1 , (13) 

where Q is any nonsingular N x N matrix and k a nonzero number. They are also invariant 
under the index reflections 



jjkl 
jykl 



tA3 
jjlk 



(14) 
(15) 



the first of these is the usual matrix transposition and the second follows from R — ► PRP, where 
P^j = Sjdj is the permutation operator. 

For special choices of Q one can obtain other discrete transformations. In general, if a is a 



permutation of the set {1, . . . , N}, then Q 
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5 % a fj\ yields a relabeling of indices Rfj —>■ R[ 



cr(k)a(l) 
a(i)a(j) ■ 



An important special case of this is given by Q*- = 5 l N _j +1 which, when followed by transposition, 
yields reflection across the secondary diagonal. 

In the present paper we restrict the i?-matrix to be an upper triangular matrix, and then 
it is natural to use only upper triangular transformation matrices Q, which we will call U from 
now on. Under such tranformations the upper triangular nature of R is preserved. Furthermore, 
one easily finds that the diagonal blocks transform as 



pfcm jri' jyk' mTj—lk 
n im ~^ u i a i'm u k' 



(16) 



Thus, in order to fix the remaining rotational freedom we will just need to impose conditions on 
the diagonal blocks. 

Normally one would require the canonical form of an upper triangular matrix be the Jordan 
canonical form, but this may require transformation matrices that are not upper triangular. For 
this reason we must relax the definition of what is a canonical form. It turns out that using just 
upper triangular transformation matrices one can bring all other upper triangular matrices into 
one of the following five 'semi-canonical' forms 





( a 


■ A 




(a 


b 


. \ 




fa 




Ci = 




b . 


, c 2 = 




a 




,c 3 = 




c . 




\ ■ 






K ■ 




c / 




\ ■ 


■ a ) 




(a 


■ A 




fa 


b 


. 








c 4 = 


\ ■ 


b c 
■ b J 


, c 5 = 


K ■ 


a 


b 

a j 









(17) 



This is the basis of our classification scheme. 
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4 How the equations were solved 

4.1 Breakdown into smaller sets 

As usual we must use the available symmetries to divide the problem into several smaller ones. In 
particular we want to fix the continuous symmetries related to upper triangular transformation 
matrices. The detailed breakdown is given in the Appendix, here we give just the general idea. 

We start with the top diagonal block and transform it into the semi-canonical Jordan form 
(|l~7|). We may therefore assume that the upper block is of the type Cj in (|l7|). We work through 
the cases in opposite order, C5 first. Using a reflection across the antidiagonal when necessary 
we may assume that the lowest diagonal block cannot be reflected up and transformed into 
anything that has already been analyzed. Using the notation that ~ C n is anything that can 
be transformed to C n using upper triangular matrices, we can write the first division by giving 
the upper and lower blocks: 

1. up C5, down anything 

2. up C4, down anything but ~ C5, 

3. up C3, down ~ C3, or ~ C4, or ~ C\, 

4. up C2, down ~ C4, or ~ Ci, 

5. up Ci, down ~ C\. 

In the third and fourth cases note that an antidiagonal reflection takes C4 into C2. 

This the starting point. Further complications arise when c = a in C3 or a = b = c in C\, 
because then they commute with an upper triangular transformation matrix having units on the 
diagonal. In those cases one can use the remaining freedom to operate on the other diagonal 
blocks. In the appendix we have given the more detailed breakdown. 

4.2 Solving the equations 

For each subcase we still have to solve the 729 equations. Some of these equations vanish 
automatically because of the upper triangular ansatz. It turns out that there are also many 
simple equations that factor to linear factors. Since the present equation solver routines cannot 
handle well very large set of equations, we decided to start the solution process by interactively 
splitting the process into branches and subbranches based on the simple factorizable equations. 
[All algebraic work was done using REDUCE pi.] The algorithm was the following: 

1. Initialize, in particular record the expression that cannot vanish (Ex: A^O). 

2. Are there simple factorizable equations? 
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Yes: Choose a simple equation (Ex. B(A — 1)(C — D)=0) and find its solutions. Go to 3. 
No: Go to 4 

3. Are there any solutions that do not break the nonvanishing condition? 

Yes: Add the first solution as an additional rule to the current branch and create a new 
branch (which inherits the current assigments and nonvanishing condition) for each 
additonal allowed branch, if any. With each later branch add the condition that 
the assignment of the previous branch cannot hold. Go to 2. (Ex: {{{B => 
0}, A}, {{A => 1}, AB}, {{C => D}, AB(A - 1)}}) 

No: Go to 5 

4. Are there equations you want to solve by hand? 

Yes: Do it, but if you have to assume something, create also a branch where the assumption 
does not hold. Go to 3. 

No: Solve the remaining equations using 'groesolve' || and output the allowed solutions. 
Go to 5. 

5. Is this the last branch? 
Yes: End. 

No: Take the next branch. Go to 2, 

A log of the interactive solution process was saved using the unix 'tee' utility, and checked 
later. Altogether the printout extended well over 1000 pages. In the end we had 552 solutions, 
but naturally most of them turned out to be of the simple type of Sec. [2|, or subcases of other 
solutions. 

5 The nontrivial solutions 

As expected, the problem of classifying the solutions is almost as time consuming as finding 
them. We must characterize the solutions by something that is invariant under the possible 
transformations. The method we finally chose was to classify first by the number of different 
eigenvalues of the i?-matrix. In counting the different eigenvalues one must be careful, because 
some solutions only seem to have the required numbed of eigenvalues. [An example is provided 
by diag(l, £i, £f , £f> £2> £l£i> £i> &)> £f = 1 ; which seems to have six different eigenvalues. 
Since there are only three different cubic roots of unity there are actually at most three different 
eigenvalues, arranged in several ways.] 
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First a word about notation. The zeroes are represented by dots, for better readability, £ 
is a cubic root of unity in general, r/ is a cubic root ^ 1, (i.e. £ 3 = 1, rj 2 = —1 — rj), e = ±1, 
other greek symbols are also roots of some polynomial equation given in the text. Symbols 
with different subscripts have independent values. Small latin letters x,y, z,p,q, k, . . . are free 
parameters. Capital letters have special properties, given in text. 

We would like to repeat again that the simple solutions of Sec. |2| are not included. Also, 
solutions related by the discussed transformations are not mentioned separately, neither are 
those obtained from previous solutions by restricting parameters. 

5.1 9 different eigenvalues 

First we have a solution with the full number of 9 different diagonal elements. 





( 1 . . 








• p ■ 


1 — qrj 2 






. . p 2 


(1 — qr])x 


(1 - qr/)(l - qrj 2 ) 






—i 2 
qp r\ 




R9.1 — 




qrj 

qp 


—q(l — qr] 2 )/x 

—r] 2 q(l — qrj) 








2 —2 

q p rj 








2—1 2 

q p if 




V ■ • • 




Q 2 



J 

Here x is a free parameter that can be scaled to 1, since it must be nonzero. The q = rjp 2 
subcase of this solution was presented in ||. 

5.2 8 different eigenvalues 

The first solution is the multiparameter version of the standard GL(3) solution, as generalized 
in [10 1 and [|ll|] to the following that has maximum of 4 free parameters, p, c, q, k. 



( 1 



P 



qp 



A 

. qk~ l 



qc 



\ 



B J 



where A, B € {1, — q}. 
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If one imposes further relations among the diagonal entries there can be other nonzero matrix 
elements elsewhere. The first such solution is 



/ 1 . 

P 



R 



■8.2 



V 



p 



X 



p V 



— 2 2 



2 —1 

-xp jX 



1-H 



1 -A* 



2 —1 



\ 



/i 3 / 



where \i is a root of (/i 3 — l)(/t/ 2 + 1) = 0. In fact i?8.2 has 8 different eigenvalues when /x 2 = — 1 
and 7 when // 3 = 1, but we keep these two branches together. If x = we get a subcase of Rs.i- 
Other solutions of this type will appear later. 
The following is a rather asymmetric solution: 





( 1 ' • 


X 


■ ■ ■ \ 




• P 

. p^ 1 


l-£ . . 












#8.3 — 




i 

. . P - 2 e 










P e . ■ 
■ p 2 . 




V • • • 




■ ■ u 



5.3 7 different eigenvalues 

The following solution is related to Rg.i- if in the diagonal of i?g.i the parameter p is restricted 
to be a cubic root of unity we can have more nonzero matrix elements: 



/ 1 . 

e 



Rl. 



xy 



1 — qrj 



(1 — qrj)x 



2^2 



QV 



xy£ 2 (l-qrj) . \ 

y£{l - qrj 2 ) 



(1 — qi])(l — qrj 1 



-yqrj 



-q(l-qr] 2 )/x 



-r] 2 q(l - qrj) 



q 2 C 2 V 2 



d 2 J 
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Here we have two scalable parameters x and y, and they can both be scaled to 1, because x ^ 
and y = leads just to a special case of Rg.i- 

A restriction on the diagonal elements of i?8.i allows 



/ 1 



R 



7.2 



P 



P^2 



V 

where A £ {1, — q}. 

5.4 6 different eigenvalues 

The following solution is related to Rg.i and Rj,\ 
( 1 • 



P V2 



A J 



R 



6.1 



V 



p 



1 



(l-q 2 )x . 



2 —i 
q p 



P 



l(l-q 2 )/x 



(1 " I 2 ) 



4 —2 



g 2 p 1 



1 / 



The p = —q subcase was found in [12]. 

The next solution is related to Rs.2'- if we have /i = 1 in i?s.2 there is other freedom available: 



R 



6.2 



V 



Here x and y can be scaled to 1, 





X 


■ y ■ \ 


• p 

. p^ 1 








p^ 1 






e 

. p~ 2 e 








p . 






. p 2 e . 
. . e J 



if they are nonzero. 
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5.5 5 different eigenvalues 

There were no new nontrivial solutions with 5 different eigenvalues. 



5.6 4 different eigenvalues 

First we have another relative of Rs.i 

( 1 



Ra. 



V 



/i 



1 — ei 



—re 



I — ei 



-i[i 2 ye 



1 — ei 



ifie 



—ie J 



where /U 4 = 1. 

Up to now the diagonal blocks have been of type C\. The next two represent a pattern that 
we will meet several times later on, the diagonal and (1,3) blocks are of type C3 (more or less) 
and in addition there are some nonzero entries on the fourth off-diagonal. 

/ 1 . x 
P ■ 



Ra.: 



Q 



P 



-1 



y 



p 



c . a \ 

. p(x + c-py -q) . 

c 



P 



X 



1 / 



/ 1 . X 

p ■ 
. 1 



Ra.; 



x-y 



p^ 1 . z 
e 

p~ l 



—x 



—xy \ 



—zp 



y-x 



P 



y 

1 J 
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( 1 • 




X 


y ■ 


. \ 






1 




■ yei 


y 








ei . 






Raa — 






• £2 • 

• ei 














p ■ 


















V • • 








p J 



5.7 3 different eigenvalues 

The next solution reduces to Rs.2 if y = 0, c.f. also -R7.1 





( 1 . . 






• -*6 2 6 2 






• 6 • 


1-6 






-y6 




• • & 






1-6 








6 2 6 






y 


R3.1 — 






66 2 


1-6 












66 












e 2 6 2 






V ■ ■ ■ 








1 





/ 1 . . 




■ y 




• £ • 




■ ■ p 




. . e 










e 2 • • 




-R3.2 — 




• • t 










s ■ 








. e ■ 




v • • • 




. . 1 
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( 1 • 




X 




^(1-6) 




• 6 


& 


■ V 












■ ■ 






#3.3 — 






• 6 • 

• • 6*6 




-1 + 6 












6 2 ? 2 




V • • 
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The following reduces to i?§.3 if y = 





( 1 • • 


X 


. . . \ 




• 6 • 


1-66 • • 






• • 6 2 


y ■ 








6 • • 




#3.4 — 




l . 

p 2 










el • • 








• 6 2 • 




V • ■ ■ 




• • 66 / 



5.8 2 different eigenvalues 

Now most of the solution have the previously mentioned structure: the diagonal and (1,3) blocks 
are like C3 and in addition there are some nonzero entries on the fourth off-diagonal. We list 
these solutions in decreasing number of nonzero elements on the fourth diagonal. 



/ 1 . x 
ei • 



#2.1 — 



V 



e\ . xe\ + ke\7}(e2 — 1) 
ei 



aei + keiUe 2 - 1) 



2^-1 



k 2 p 



1 / 
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1 . y — x 

e 

. . 1 


• P 

—x 


-y + x 


-ye 


\{-2y 2 +x 2 ) 
-(y + x) 




e . ye 
. -1 

e 




x 








1 


€ 


y + x 
1 



R 



2.3 




/ 1 



R 



2.4 



X 



P 

x — k 



e . ke 
. -1 

e 



ae . q 
a 

ae 



1 . x 

e 

. . 1 / 



/ 1 • 



#2.5 = 



ei • 



P 



ei . xei 
ei 



— x . — xy \ 
—xe\ 

-y 



1 7 
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For the following two the diagonal blocks are of type C2. If the parameter b were zero one 
could make an index shuffle so that they would fit into the previously mentioned band form. 





fix 


2 

— X — X 






(I 


\ 




1 


— X 






-b 






■ ■ ei 




— Xt\ 












1 X 






b 




R2.6 — 




1 


£1 
















€\ X€\ 














£1 








V • • • 








£2 


/ 





tlx. 

. 1 . 

■ ■ ei 


—x x(x — 2c) 
x - 2c 

— cei 


. -6x \ 
-6 






1 -x + 2c 




R2.7 — 




1 








ei 






V • • • 




ei cei 
ei 

• • £2 / 



The following are rather sparse, but have more freedom on the diagonal 





( 1 • • 

• ei • 

• • ei 


x a 


. b c \ 






ei . 




R2.8 — 




• £2 • 








• £2 






V • • • 




£1 • 

• £2 • 

• • £2 / 
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( 1 • • 

• ei 




. . a \ 






£3 • 




-R2.9 — 




. e 4 . 














V • • • 




ei . 

• • e 5 J 



5.9 One eigenvalue 

The next four solutions have basically the C3 structure, with some extra elements. 





( l 


1 


. . -y + x 


—x . — xy 
—x q 

-y 








1 . 


X 


y-x -q 


R\.i — 




















1 






V ■ • ■ 




1 • y 

1 

1 





( 1 


X 

1 


• • p-q 


p ■ 


a 
k 

P 


\ 








1 . X 




-k 




R1.2 — 






. 1 

. . 1 
















1 . 


X 












. 1 








V • 








1 


/ 
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Rl.3 — 


/ 1 . X 
. 1 
. . 1 




—x . —xy \ 

-y -p 
-y 






1 • y 

. 1 . 
. . 1 


p 

q 


RlA = 


f 1 


V • 

. X 
1 . 
. 1 






i • y 
l 

• • i ) 

—x 

. -2y + x 


-xy ^ 
P 

-y 




1 . 2y - x 
. . 1 




V • • • 




1 • y 
1 

1 J 



The next one is like C2. The labeling change 2^3 would make it almost C3 but would not 
keep the solution upper triangular. 





fix. 




—x —x 


—k + xq 


p k 


a 


\ 




. 1 . 


—x 






b 






. . 1 




—x 




P 








1 X 




■ p-q 


-b 




-R1.5 — 




1 


1 
















1 x 














1 








V • • • 








1 


/ 



The next three solutions have mixed diagonal blocks. It could be possible to transform the 
upper block to C2 or C3 but we chose a different and apparently simpler form defined by the 
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requirement that the above mentioned relabeling 2^3 keeps the matrix upper triangular. 





/ 1 x p 
. 1 . 
. . 1 


—x 


-xy 

-y 


— 2px — k 
—x 


— p k 
■ ~P 


-pq 

-q 






1 


y 


p 






-Rl.6 — 






l 
















1 








V • • • 




1 X 

1 


q 
l 





( 1 x p 

. 1 . 
. . 1 


— X 


-xy 

-y 


c 

-q+p 
—x 


a 


— c + x(—q + a) b \ 
a 

a 






1 


y 


q 






Rl.7 = 






l 
















l 








V • • • 




1 


x p 
1 

1 / 





f 1 x q 


—x 


-x 2 


-k 


— q k 


px 


\ 




. 1 . 
. . 1 




—x 


q 

—x 




p 




-Rl.8 — 




1 


X 

1 


1 


■ -q 


-p 














1 X 

1 








V • • • 










1 


/ 



Finally we have two solutions whose diagonal blocks are of type C5. 





/ll. 


-1 -1 X 


y 


y-x 


p 


\ 




. 1 1 


. -1 -1 




y 


-z + y 






. . 1 


. -1 






y 








1 1 


-1 


-1 


z 




#1.9 — 




1 1 




-1 


-1 








1 






-1 










1 


1 














1 


1 






V • • • 








1 


/ 
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In the next one the parametrization is rather complicated and it would be interesting to 
understand its origin, [x = 1 yields a symmetric solution.] 



/ 1 2 . 
1 2 
1 



Ri 



10 



-2 -Ax 8x 
. -2x 4 - 6x - 2x 2 
2 -Ax 



1 2x -2x(l - x) 
1 2x 
1 



4 . -4x(7 + x) \ 

2x(l + x) 6x(2 + x)(-l + x) 
-Ax(l - 2x) 



-2 -2x(l + x) -2x(l + 3x)(-2 + x) 
-2x 4x(l - 2x) 

2 -Ax 



1 -2 + 4x 4(l-2x)(l-x) 
1 -2 + 4x 

1 



6 Discussion 



In this paper we have given all upper triangular nonsingular solutions to the constant quantum 
Yang-Baxter equation (Q), modulo upper triangular transformations ( |i~3"l) and reflections (14,15) 
and minus simple solutions of Sec. ||. 

There is not much that we can say now about the solutions, here are just some random 
observations, (i) Some of solutions (-R9.1, R%.\ — Rs.3, R7.1, R&.i-, R&.2, R4.1 and R3.1 — R3.4) 
satisfy also the weight condition R^j = 0, if k+l 7^ i+j (mod 3). Several others satisfy the same 
condition (mod 2), e.g. R2.1 — i?2.5- (h) The upper triangular property is of course sensitive to 
index relabeling, for example 2^3 often breaks it. Sometimes one can find an upper triangular 
transformation after which the solutions stays upper triangular even after 2 <-> 3 exchange. 
This happens for the solutions Rs.3, R4.2, Raa, R3.4, R2.4, R2.5, R2.9, R1A1 R1.6, R1.7 an d R1.8, 
and they are presented in this special form, (iii) Some solutions show interesting parametric 
relations, in particular Rg.%, R7.1, R6.1 and -Ri.io- 

Now that we have these constant solutions several natural questions come up for further 
study, for example, can one add a spectral parameter, can one contruct a corresponding universal 
R-matrix? The corresponding algebraic structures also need investigation. 



Note for the readers of the preprint version: The bibliography may contain omissions, if you know 
additional references where explicit three-state upper triangular solutions have been presented, 
please let me know by e-mail to hietarin@utu.fi 
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Appendix 



In this appendix we give the breakdown of the original problem to smaller subsets. The primary 
classification proceeds by the upper triangular blocks on the diagonal. The "number of solutions" 
given below is the number before any trivial solutions or subcases were eliminated. 

A Upper block of type C5 

We use scaling freedom to put all nonzero entires in C5 to 1. The solution procedure produced 
8 solutions, which were immediately reduced to 3 basic solutions. 

B Upper block of type C4 

Let us scale so that b = c = 1 in C4. In principle we could exclude lower blocks of type C5 from 
the very beginning, but no such solutions were found anyway. 

1. a arbitrary, C5 7^ or f\ 7^ 0: 3 solutions, all have a = 1. 

2. a ^ 1, C5 = /V = 0: 22 solutions 

3. a = 1, C5 = /V = 0: 1) 67 7^ 0, 4 solutions. 2) 67 = 0, a§ = 1, 9 solutions. 3) 67 = 0, ag = 0: 
35 solutions. 

C Upper block of type C3 

We can scale a = b = 1 

C.l c 5 ^0or/ 7 ^0 
no solutions 

C.2 c 5 = / 7 = 0,c^l 
33 solutions 

C.3 c 5 = / 7 = 0, c = 1 

Now C3 commutes with any upper triangular matrix with units on the diagonal. Using this we 
can put the lower block also into a semicanonical form. Solving some of the equations reveals 
that the diagonal entries of the lower block are also all = 1. If the lower block is of type C5 a 
reflection takes the system to one studied before, same holds if the lower block is of type Ci. If 
the lower block is of type C4, we get one solution. 
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What remains is a lower block of type C3 with a = c, but possibly with 6 = 0. This again 
commutes with UT transformation matrix with units on the diagonal, and thus we will use this 
rotational freedom to put the center block in a semicanonical form. 

1. Center block of type 5: no solutions. 

2. Center block of type 2: 3 solutions. 

3. Center block of type 3: 42 solutions. 

D Upper block of type C2 

Now lower block must be of type C4 or C\, all others can be reflected to cases studied before. 

1. /V 7^ 0: No solutions. 

2. C5 / or f 7 = 0: 2 solutions. 

3. C3 7^ 1, C5 = fr = 0: 11 solutions 

4. C3 = 1, C5 = ff = 0: In this case one quickly finds that the lower block must be 



Here we need to analyse only the case h = g = due to the following: 1) If k = I, h 7^ 
we can use using upper triangular transformations to put g = 0. Then lower block is of 
type C2 which reflects to an upper block of type C4 done earlier. 2) If k = I, h = 0, g 7^ 
lower block is of type C3, done earlier. 3) If k 7^ I, fe/Owe can transform lower block to 
type C2 without changing the upper block. 4) In the remaining cases the lower block is 
diagonal, there were 15 solutions of this type. 

E Upper block is nonunit diagonal, C5 7^ 0, fj ^ 

Let us normalize the diagonal elements as (1,6, c). Here we may also assume that the lower 
block is diagonalizable by an upper triangular matrix, but maybe not yet in diagonal form. In 
most cases solving the first few equations yields directly a diagonal lower block, exceptions to 
this are discussed separately. 

In general one finds that if C5 / or fa / then C3 = 6|, this implies in particular that in 
this section we must have 62 7^ 1. After that one also finds quickly that the lower block must be 



/ k 



h 



k 



v • 



1 J 



( h 



9\ 



k 



(18) 



V • 
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1. b 3 ^ 1,6 2 ^ 1: 3 solutions. 

2. 6 = — 1: As mentioned above we must now have C3 = 1. If now I = h in ( |l8|) the system 
can be reflected to one of the cases studied before. If I 7^ h one can use the transformation 

(1 ■ y\ 
. 1 . , 

V • • 1 / 

which commutes with diag(l, —1, 1), to put g = 0. One finds 4 solutions. 

3. b 3 = 1, b 7^ 1: 5 solutions. 

F Upper block is nonunit diagonal, C5 7^ 0, fj = 

The same comments as above holds in this case and we just give the results. 

1. b 3 ^ l,b 2 ^ 1: 4 solut 10ns. 

2. b = —1: 4 solutions. 

3. b 3 = 1, b 7^ 1: 6 solutions. 

G Upper block is nonunit diagonal, C5 = 0, 7^ 



As before, the lower block must be as in (|18|). If Z = h we can reflect the system to case F or C, 
else we eliminate g and reflect to case F. 

H Upper block is nonunit diagonal, C5 = 0, f-j = 
H.l The diagonal elements (1,6, c) are all different 

76 solutions 

H.2 c= 1,6^1 

Again one finds that the lower block must be of type (|l^) and as in ^.2 one can show that the 
lower block must in fact be diagonal. If its diagonal elements are all different we reflect it to 



H.l . Thus we get subcases according to the diagonal entries of the lower block: 

1. Lower block = diag(l, k, I): 21 solutions 

2. Lower block = diag(k, k,l) or diag(k,l,l): 37 solutions. 

3. Lower block = diag(l, I, I): 19 solutions. 
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H.3 6=1, c^l 

Now we may assume that the lower block is transformable to a diagonal form, which has at 
most two different elements and is not (h,l,h), all other cases can be reflected to ones studied 
before. Then one finds that the lower block must be 

(h g .) 

. k . . (19) 
{. . l) 

Since the upper block stays invariant under transformations 

(i y ■ \ 
. l . , 

V ■ • i/ 

we can transform the case h ^ k to a diagonal form while the case h = I reflects to C4. 

1. Lower block is diag(h, h, I): 36 solutions. 

2. Lower block is diag(h, I, I): 18 solutions. 

3. Lower block is diag(l, I, I): 36 solutions. 

H.4 b = c^l 

Again one can trasform the lower block to something that reflects to one of the cases before, 
except for two cases: 

1. Lower block is diag(h, h, I): 17 solutions. 

2. Lower block is diag(l, I, I): 40 solutions 

I Upper and lower blocks are unit matrices 

In this case we have the full trasformation freedom left and we use it to bring the center block 
into the semicanonical form (recall C2 reflects to C4) 

1. Center block type C5: no solutions. 

2. Center block type C2: 2 solutions. 

3. Center block type C3: 13 solutions. 

4. Center block type C±: Let us call the diagonal elements of the center block as (d, f,g). 1) 
d, f,g all different, 4 solutions. 2) g = d, f ^ d, 6 solutions. 3) d = f, f 7^ g, 9 solutions. 
4) Center block unit matrix, 9 solutions. 
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